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The effects of varying the strengths of tensor and spin-orbit interactions on Ml and 
E2 rates in 12 C: Comparison of results in AiV = and = + 2hu> spaces 

M.S. Fayache, 1 ' 2 , Y.Y. Sharon 2 ! and L. Zamick, 2 
(1) Departement de Physique, Faculte des Sciences de Tunis, Tunis 1060, Tunisia 
(2) Department of Physics, Rutgers University, New Brunswick, NJ 08903 USA 

The energies and transition rates to J — 1 + T = 1 and J = 2 + T = 0,1 states in 12 C with matrix 
elements fitted to realistic G matrix elements obtained in non-relativistic approaches are studied. 
Then the effects of varying the strengths of the two-body tensor and spin-orbit interactions are also 
considered. The calculations are done in both a small space (Op) and a large space (Op + 2hu>). 
In the small space the B(M1) from ground to the J = 1 + T = 1 is enhanced and gets closer to 
experiment if the strength of the spin-orbit interaction is increased and/or if that of the tensor 
interaction is made weaker. In a large space the spin B(M1) gets reduced by almost a factor of two. 
A 'self-weakening' mechanism for the tensor interaction which succeeded in explaining anomalies in 
other nuclei does not seem to work for this case. 



I. EXPERIMENT AND PREVIOUS STUDIES 

Of necessity, shell model calculations must be carried 
in limited shell model spaces. It is therefore of interest 
to try to understand the impact of such restriction on 
the results that are obtained. In this paper we study the 
effects of enlarging the shell model space on Ml and E2 
excitation rate and we will also study the effects of vary- 
ing strengths of the spin-orbit and tensor interactions. 

There are two well studied J = 1 + states in 12 C [Q. 
These are the T = states at 12.71 MeV and the T = 1 
state at 15.11 MeV. The B(M1) rates to these states from 
ground are respectively 0.040(3) [i% and 2.63(8) fi 2 N . The 
dominance of the isovector over the isoscalar B(M1) an 
easily be understood by taking the ratio of the squared 
isovector and isoscalar spin g factors 
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Nuclear structure effect make the ratio even larger, 
ss 336 0|. Because the transition to the dominantly 
J = 1+ T = state is so weak it is strongly affected 
by isospin impurities, due mainly to the Coulomb inter- 
action ||-H ■ This has previously been considered by Sato 
and Zamick Q , where it was shown that the surprisingly 
large Coulomb matrix element needed by Adelberger et 
al. |s| to induce the isospin mixing could be explained 
by taking into account nuclear binding energy effects. In 
this work, however, we will focus of the J = 1 + T = 1 
states and so we use the experimental value of 2.63(8) fi 2 N 
as a point of comparison for our calculations. 

There is some interest in correlating B(M1) transi- 
tions both orbital @-§ and spin with B(E2)'s 
i.e. with quadrupole collectivity. We note that there is 
good evidence that 12 C is strongly deformed (/? w —0.6). 
The magnitude of (3 comes from the experimental B(E2) 
0i -> 2i which equals - 40 e 2 fm 4 

Using phenomcnological interactions Cohen and Ku- 
rath |14|] had no difficulty in fitting this data. They 
obtained the states at 12.43 MeV and 15.23 MeV with 



B(M1) values of 0.0148 fx 2 N (Coulomb interaction ne- 
glected) and 2.510 fi 2 ^ respectively. However, calcula- 
tions by Zheng and Zamick |l5)| , using G matrices calcu- 
lated by non-relativistic methods from realistic interac- 
tions, gave much too small a B(M1) to the J = 1 + T = 1 
state at 15.11 MeV. It was noted that increasing the spin- 
orbit strength and/or weakening the tensor interaction 
strength could cause the B(M1) to be increased. A justi- 
fication for increasing the spin-orbit strength comes from 
the Dirac phenomenology |l(J. Indeed, Zheng, Zamick 
and Miither |l7j] showed that with G matrices derived 
from the relativistic Bonn A interaction p8| , the value 
of B(M1) + T = 0^1+T=l increased with de- 
creasing Dirac mass. For example, with to* (DIRAC) = 
938.9, 729.1 and 630.0 MeV/c 2 respectively, the values of 
B(M1) were 0.69, 1.13 and 1.57 n 2 N respectively. Zheng 
and Zamick had introduced the (x, y) interaction |15j 



V = V c + xV s . 



yv t 



(1) 



where s.o. = spin orbit, t = tensor and c = everything 
else, especially a spin-dependent central interaction. For 
x = 1, y = 1 this interaction gave a good fit to the non- 
relativistic Bonn A matrix elements. On could simulate 
the effects of the Dirac phenomenology very well by mak- 
ing x larger than one. 

One can also get a larger B(M1) by making y smaller. 
The justification for this is less clear. There are the uni- 
versal scaling ideas of G. E. Brown and M. Rho [[l9| where 
inside the nucleus the masses of all mesons except for the 
pion decrease relative to their free space values. In this 
model not only does the spin-orbit interaction become 
stronger (as in Ref. |l(|), but also the tensor interaction 
becomes weaker. Another possibility is the that the effect 
could be due to higher shell admixtures, without having 
to invoke non-nuclear degrees of freedom. Concerning 
the latter, Fayache, Zheng and Zamick |2rj| noted that 
higher shell admixtures, i.e. 2Tiuj and higher excitations, 
could for many phenomena make the tensor interaction 
look weaker in the valence space. They called this the 
'self weakening mechanism'. They applied this to the 
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problems of the quadrupole moment of the ground state 
of 6 Li, the near-vanishing of the Gamow- Teller matrix 
element for the A = 14 system ( 14 [J = T = 1} 
-> 14 N [J = 1 T = 0]) and the energy splitting of the 
J = 0^T=landJ = 07T = states in 16 0. In 
all three cases the experimental result could be obtained 
either by weakening the tensor interaction in a valence 
space calculation i.e. making y less than one, or by keep- 
ing the full strength y = 1 and performing shell model 
calculations in which 2tku and sometimes higher admix- 
tures were included. 

However, in a perturbation theory approach by Zheng, 
Zamick and Muther |[7| to the problem of B(M1) 12 C 
[J = T = -» J = 1+ T = 1] the self-weakening 
mechanism did not seem to work. The authors renor- 
malized the two body matrix elements in the Op shell and 
then used these in a Op shell matrix diagonalization. If 
only Bertsch-Kuo-Brown bubble pl] , p2|| was exchanged 
between two nucleons, then indeed the B(M1) got big- 
ger. However, this was offset by the two other diagrams 
- hole-hole and particle-particle The latter re- 

inforce each other and cause B(M1) to become small 
again The specific number for the Bonn A interaction 
were 0.639 /ify for the bare interaction, 1.177 n 2 N when 
the bubble diagram was included, 0.832 ii 2 N when the 
bubble plus hole-hole were included and 0.514 ^l 2 n when 
the bubble plus hole-hole plus particle-particle diagrams 
were included. We go from small to large to smaller to 
smaller still. 

It is not however clear if merely renormalizing the two 
body interaction is sufficient for calculating a transition 
matrix element. If one draws Goldstone diagrams for the 
Ml transition there are several diagrams not included by 
merely substituting a renormalized interaction into the 
shell model calculation. Examples of the latter include 
the diagrams where the Ml operator is inserted after the 
first interaction but before the second interaction. There- 
fore in this work we decided to do a complete valence plus 
2huo shell model diagonalization for the relevant states 
and the Ml transition rates in 12 C. In such an approach 
all the missing diagrams are implicitly generated. 

Of particular interest is whether the large space calcu- 
lation will or will not serve to justify the adjustment of 
the parameters x and y in the small space. For example, 
will and x = 1, y = 1 calculation in the large space yield 
results similar to those of an 2; = 1, y = 0.5 calculation 
in the small space (self-weakening of the tensor interac- 
tion)? Such a comparison makes sense only if we look at 
the low energy sector of the large space results. 

Our calculations for the B(Ml)'s are of additional 
interest because of the close relationship of the spin 
B(Ml) CT 's to the beta decay Gamow- Teller matrix ele- 
ments. Thus our studies of the effects of enlarging the 
shell-model space have a broader applicability. 



II. COMPARISON OF RESULTS OF B(M1) 

J = 0+ T = -> J = it T = 1 IN 12 C IN THE SMALL 
SPACE (AN = 0) AND THE LARGE SPACE 

(AN = 2) 

A. The x = 1, y = 1 results 

In the small space (see table ||) with the x = 1, 
y = 1 interaction (which corresponds roughly to G ma- 
trix elements of a realistic interaction obtained in a non- 
relativistic approach) we obtain for B(M1))| to the low- 
est J = 1 + T = 1 state a value of 0.89 /i^, much smaller 
than the experimental value of 2.63(8) fi%. Our calcu- 
lation further indicates that this is dominantly a spin- 
mode excitation, since B(M1); = 0.03 rfL and B(M1) (T = 
0.58 y? N . 

Does the discrepancy between theory and experiment 
go away when we go to the large space? Quite the con- 
trary, the B(M1) gets smaller, mainly because B(M1) CT 
decreases. Note the B(M1) CT is very closely related to 
B(GT) the Gamow- Teller transition. They will have 
identical quenching due to configuration mixing. From 
table I we see that relative to AN = 0, the quenching 
factor for B(M1) CT in a AiV = + ITao calculation is 
0.31/0.58 = 0.53 (or for the operator it is 0.73). 

This is in accord with the general consensus of what 
the quenching of Gamow- Teller operators should be. We 
have here verified this result, not in perturbation the- 
ory, but rather in an explicit matrix diagonalization. It 
appears that matrix diagonalization results are close to 
those of perturbation theory. 

A point that has been made before jl5| is worth re- 
peating. In open shell nuclei, one can get B(Ml) CT 's and 
B(GT)'s to be too small, despite the fact that we believe 
there should be an overall quenching. This is what is 
happening here when we use non-relativistic G matrixes 
derived from realistic interactions. As another example, 
in the SU(4) limit for N = Z nuclei B(GT) will vanish. 
This limit can be reached by using a spin-independent 
two-body interaction and turning off the one-body spin- 
orbit interaction, e.g. in the Elliott model where one has 
only a Q ■ Q interaction j2j|. Thus in a given calcula- 
tion in an open shell nucleus one cannot readily deduce 
where one needs quenching or enhancement of the spin, 
magnetic and/or Gamow- Teller operator, unless one has 
complete confidence in the two-body interaction and sin- 
gle particle energies that are being used. 

In table |l](a) we present the summed Ml strengths 
in the small space; in table 0(b) we have the summed 
strengths to the first 10 states in a large space calculation 
(AN = Q + 2%u) and in table f|(c) the sum to 500 states. 

Whereas for the lowest state in the small space calcula- 
tion (table |) for the x = 1, y = 1 interaction the values of 
B(M1), B(M1); and B(M1) (T were respectively 0.89, 0.03 
and 0.58 the summed small space values are 1.42, 
0.60 and 0.91 n 2 N . In particular there is some isovector 
orbital (scissors mode) strength in higher J = 1 + T = 1 
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states. 

Focusing on B(M1) CT we see that the small space sum 
of 0.91 fj? N becomes 0.70 [i 2 N when we consider the first 
10 states in a large space calculation. There is quenching 
of the low lying strength. We will compare the small and 
large space total strengths in a later section. 

B. Varying x and y 

In tables [| and [El] we study also the variation of the 
B(M1) rates with x and y, the strengths of the two- 
body spin-orbit and tensor interactions respectively. We 
consider four sets of (x,y) - (1,1), (1.5,1), (1,0.5) and 
(1.5,0.5). Increasing the spin-orbit strength from x = 1 
to x = 1.5 simulates to a large extent the use of Dirac 
phenomenology with a Dirac effective mass of 1/1.5 = 
0.67. If indeed the source of the increase in the spin- 
orbit interaction comes from the Dirac phenomenology, 
then one is justified in using this enhanced value of x 
in both the small space (AN = 0) and the large space 
(AN = + 2hw). 

For the tensor interaction the choice of the parame- 
ter y in the large space calculation depends which of the 
scenarios discussed in the previous section for weakening 
the tensor interaction are correct. If the universal scaling 
ideas of G. E. Brown and M. Rho |ll| are correct then 
one should use a weaker tensor interaction, e.g. y = 0.5 
in both the small space and large space calculation. How- 
ever in the 'self-weakening mechanism' which explained 
several phenomena for nuclei with mass numbers A = 6, 
14 and 16, the higher shell admixtures were responsible 
for making the tensor interaction appear weaker in the 
valence space. In that case if one uses y = 0.5 in the 
small spaces one should use y = 1 in the large space and 
hope that the results for the low energy properties are 
about the same in the two calculations. 

First, focusing on the small space calculation of table | 
we see that we get a marked improved in the B(M1) rate 
J = + T = 0^J=1 + T=1 when we increase the 
spin-orbit interaction and/or decrease the tensor interac- 
tion. From (x,y) = (1, 1) to (1.5,0.5) the value of B(M1) 
increase from 0.89 to 2.54 and the value of B(M1) CT in- 
crease from 0.58 to 2.57. We are getting close to the 
experimental value of 2.68 n 2 N . 

In the large space the values of B(M1) to the first 
1 + T = 1 state are however in all cases smaller than 
in the small space. The 'self- weakening mechanism' for 
the tensor interaction does not appear to work here. We 
should compare the (x,y) = (1.5,1) large space B(M1) 
result of 1.29 ^l 2 n with the (x,y) = (1.5,0.5) small space 
result of 2.54 y? N . For (x, y) — (1.5, 1) the gains from the 
small space to the large space does not increase B(M1) 
- rather it decreases it from 1.89 fi 2 N to 1.29 /i^. This 
is mainly due to a decrease of B(M1) CT from 1.71 n 2 N 
to 0.98 [1% ( B ( M1 )/ changes from 0.04 to 0.02). Why 
the 'self- weakening mechanism' appears to work in some 



nuclei but not in others is not clear and bears further 
investigation. 

It should be added that results we obtain here via ma- 
trix diagonalization are not so different from those by 
Zheng, Zamick and Miither [ fL7| using a renormalized in- 
teraction in perturbation theory. Evidently, the missing 
diagrams in which the Ml operator acts between the two 
interactions are not so important. 

We see from table [n[ which deals with summed B(M1) 
strengths, that in the small space increasing the spin- 
orbit strengths from s = ltoaf = 1.5 results in very 
significant changes in the summed isovector Ml excita- 
tion strengths. More specifically, with this change in x, 
B(M1) is increased by over 60%, B(M1) CT is increased by 
127% and B(M1) Z is decreased by 16%. The effect in the 
small space of decreasing the strength of the tensor force 
from y = 1 to y = 0.5 is similar but smaller; B(M1) is in- 
creased by 11%, B(M1) (T is increased by 24% and B(M1)/ 
is decreased by 1.6%. The combined effect of increas- 
ing the spin-orbit and decreasing the tensor strengths is 
more than the sum of the individual effects. Here B(M1) 
is doubles, B(M1) CT is more than tripled which B(M1)/ 
decreases by 23%. 

From table || we also see that in the large space 
the effect on the isovector B(M1) strength of increasing 
strength x of the spin-orbit force is qualitatively similar 
to the effect of doing so in the small space, but is smaller 
in scale by a factor of about 0.6. When x is increased 
fro 1 to 1.5 (both with y = 1), the B(M1) is increased 
by 36%, B(M1) CT is increased by 80% and B(M1) Z is re- 
duced by 10%. Very interestingly, in the large space, 
changing only the strength y of the tensor force has es- 
sentially no effect on the summed isovector Ml strength. 
Indeed, when x — 1 and y is changed for 1 to 0.5, then the 
B(M1), B(M1) CT and B(M1); all change by less than 3%. 
This result requires further study. On the other hand, 
in the large space, the effect on the isovector B(Ml)'s 
of decreasing the tensor strength y (from 1 to 0.5) af- 
ter the spin-orbit strength has already been increased for 
x = 1 to x = 1.5 is much more substantial. It results in 
a further increase of 17% in B(M1) and close to 30% in 
B(M1) CT and a further decrease of 8% in B(Ml)j. 

The middle part of table || shows for any given (x, y) 
combination what happens if we limit ourselves to the 
first ten states in the large space. Then we typically un- 
derestimate the isovector B(M1), B(M1) CT , B(M1) Z val- 
ues by 10% to 30% in comparison with the results ob- 
tained by summing over the lowest 500 states in the large 
space. In the small space there are only eight states with 
J = 1 + T = 1. For any one (x,y) combination the cor- 
responding values of B(M1) and B(M1) (T are typically 
20-40% larger than for the ten lowest states in the large 
space, but the B(Ml);'s hardly change (by less than 7% 
in all cases. 
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C. Comparison of the total isovector Ml strength in 
the small and large spaces 

By including 2Tiuj excitations in our shell model diago- 
nalization we are able to address the question of whether 
or not the isovector Ml strength which disappears from 
the low lying sector (and in particular from the lowest 
1 + T = 1 state at 15.11 MeV) reappears in the high en- 
ergy sector, i.e. in the 2huj region. To answer this we 
merely have to look at the results in table ||. 

For the case x = 1, y = 1 the total strength in the 
0p+2?ia; space is larger than in the small space. For the 
other three cases (x,y) = (1.5, 1), (1, 0.5) and (1.5,0.5) 
the opposite is true - the total isovector magnetic dipole 
strength is smaller when Ihio excitations are included. 

For x = l,i/ = 1 the summed B(M1) strength in the 
small space is 1.42 fj,% whereas in the large space it is 
1.61 In this case, where there is quenching in the 
low-lying sector this is more than compensated for by 
the strengths in the 2hu> region. 

For the other three cases, the reverse is true. For ex- 
ample in the case x = 1.5, y = 0.5 the summed strength 
in the Op space is 2.58 fj^, whereas in the 0p+2hco space 
it is 2.56 n 2 N . If we look at B(M1) CT the drop in the sum 
is larger, from 2.82 ^ 2 N to 2.28 a decrease of 19%. 
The orbital summed strength B(M1); on the other hand 
increases for this case from 0.46 fi 2 N to 0.57 fi 2 N . Indeed 
the orbital summed strength increases for all four (x, y) 
combinations considered here. 

From the above we see that the Ml strength is redis- 
tributed but not conserved. There is definately a quench- 
ing in the low energy sector but in three out of four cases 
considered here the missing strength does not reappear 
in toto in the high energy sector. 



III. COMPARISON OF RESULTS OF B(E2) IN 
12 C IN THE SMALL SPACE (AN = 0) AND THE 
LARGE SPACE (AN = 2) 



As a counterpoint to the high lying Ml transitions we 
wish to now study what happens to low lying states, e.g. 
the 2^~ state, when we increase the model space and when 
x and y are varied. In the small space (Op) we use ef- 
fective charges e p = 1.5, e n = 0.5. In the large space 
(0p+2?iu;) we use bare charges e p = 1, e„ = 0. The 
AiV = 2fwj admixtures should in principle at least jus- 
tify the use of effective charges in the small space. The 
results of the B(E2)'s to individual states are presented 
in table III and the summed strengths in table IV. 



A. The x 



1 results 



1. Isoscalar B(E2) -J = 0+T = 0^J = 2+T = 

In the small space there are eight J = 2 + T = states. 
In that space (with the effective charges e p — 1.5, e n = 
0.5), 95% (80.6 e 2 fm 4 ) of the total calculated B(E2) 
strength of 82.7 e 2 fm 4 is concentrated in the transition 
to the lowest J = 2+ T = state at 3.80 MeV. In the 
larger space, the lowest J = 2 + T = state now at 4.79 
MeV is still dominant with a B(E2) of 23.2 e 2 fm 4 and 
accounts for 55% of the total strength of 42.5 e 2 fm 4 . In 
addition, another relatively strong J = 2 + T = state 
emerges among the ten lowest state, its energy is 4.40 
MeV and its B(E2) is 6.12 e 2 fm 4 . The B(E2) isoscalarT 
sum, over the first 500 states in the large space, is slightly 
less than one half of the total B(E2)| sum over all the 
states in the small space calculation. The sum over the 
first ten states in the large space (up to about 45 MeV 
of excitation) adds up to 31.9 e 2 fm 4 or about 3/4 of the 
sum over the first 500 states. 

The dominance of one single state suggest collectivity. 
In N — Z nuclei the isoscalar B(E2) is proportional to 
(e p + e„) 2 . Usually the values of the effective charges 
e p , e„ are introduced in AN = calculations to repro- 
duce the effects of neglecting higher shells. The com- 
parison of the large space and small space results sug- 
gests that e p + e n should be more like 1.4 instead of 
1.5 + 0.5 = 2. On this point it was noted by Abbas 
and Zamick H| that the isoscalar effective charges is 
less than two for light nuclei and only reaches the value 
of two asymptotically for large A. 



2. Isovector B(E2) - J = + T = 0^J = 2+T = 1 

In general this isovector B(E2) sum is smaller than the 
corresponding isoscalar one. In the small space, with a 
total of seven J = 2 + T = 1 states the isovector B(E2)j 
sum is dominated by the second J = 2 + T — 1 state 
(here at 15.78 MeV), which contributes 2.61 e 2 fm 4 out 
of a total value of 3.39 e 2 fm 4 for the B(E2) sum. In the 
large space calculation the largest contribution to the 
isovector B(E2) sum from among the few lowest- lying 
state again comes from the second excited J = 2 + T = 1 
state now at 20.08 MeV above the ground state with a 
B(E2) of 1.56 e 2 fm 4 . The sum over the B(E2) T 's to 
the first ten J = 2 T = 1 states (which corresponds to 
an excitation energy of up to 53 MeV) is only 2.91 e 2 
fm 4 , but when the 500 lowest states are taken into ac- 
count the B(E2) sum increases about seven-fold to 21.28 
e 2 fm 4 . It would thus seem that there are substantial 
isovector B(E2) contribution from very high lying states. 
The isovector B(E2)'s for N — Z nuclei are proportional 
to (e p — e„) 2 ; reductions in both e p and e n will not neces- 
sarily reduce this quantity, unless such reduction are car- 
ried out in such a way that the difference e p — e n is also 
reduced. There is evidence that indeed e p — e n should 
be less than unity, from theoretical calculations [^5|j26| 
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and from phenomenological fits [27 2 
is 0.7. 



A popular value 



B. Varying x and y 

1. Isoscalar B(E2) - J = + T = 0^J = 2+T = 



We see from table II] the dominant B(E2)'s, the 
isoscalar ones, depend only weakly on the particular 
(x, y) combination that is used. For each of the three 
main blocks in table |l| (small space, large space, the low- 
est states in the large space) the variation of the isoscalar 
B(E2) among the different (x,y) is less than 9%. Typi- 
cally , with x = 1 decreasing the tensor force from y = 1 
to y = 0.5 has a negligible decreasing effect (less than 
1%). Increasing the spin-orbit from x = 1 to x = 1.5 
both with y = 1 decreasing the isoscalar B(E2) by 3-7%, 
and simultaneously decreasing the tensor and increasing 
the spin-orbit decreases the isoscalar B(E2) by a total of 
5 to 9%. 

The isoscalar B(E2) does depend critically, however on 
the space in which the calculation is carried out and other 
effective charges used. The approximate B(E2) values in 
e 2 fm 4 are 90- for the small space (with e p — 1.5, e n = 
0.5), 30 for the first ten 2 + states in the large spaces and 
40 for the large space. We thus see that the isoscalar 
B(E2) sum decreases as we calculate it in larger spaces 
(with the actual physical charges e p = 1, e n = 0). 

Several conclusions can be drawn from the results 
which were outlined in the previous paragraph. The 
isoscalar B(E2)'s for N — Z nuclei are proportional to 
(e p + e n ) 2 . We use effective charges e p = 1.5 and e„ = 0.5 
in the small space to compensate for the non-inclusion 
of higher shells. These charges lead to B(E2)'s that 
are twice as large in the AiV = calculation as in the 
AiV = 2 calculation. If we use the large (AN = 2) space 
calculation as our standard, it would seem the effective 
charges e p + e n = 2 used in the small space are too large 
and that this sum should be smaller by a factor of about 



1.4 (i.e. e p 



1.4) at least for C and possibly more 



general in the Op shell. 

The use of effective charges that are too large in the 
small space (AN — 0) calculations leads to exaggerated 
B(E2)'s and hence to exaggerated deformations. Such 
a calculation makes 1 C appear more deformed than it 
really is. 

In the small space at least 96% of the dominant 
isoscalar B(E2) sum is due to one state, the lowest 
J = 2 + T = state which is calculated to have an ex- 
citation energy of about 3.8-4.1 MeV (for the different 
(x,y) combinations). 

In the large space that state moves to 4.6-4.9 MeV (for 
the different (x, y) combinations) and its isoscalar B(E2) 
value increases but only 10-12%; however, now for all the 
(x, y) combinations there is another strong state (with 
strength 20-30% that of the first excited state), which 



lies high at 42-44 MeV, and is the eight or ninth excited 
J = 2+ T = state. Together, these two strong states 
supply about 70% of isoscalar B(E2)j strength summed 
over the first 500 J = 2 + T = states in the large space. 



2. Isovector B(E2) - J = + T = 0^J = 2+T=1 

Here there is a marked difference in the isovector B(E2) 
values between small and large spaces. The small space 
isovector values in e 2 fm 4 range from about 3-5 (depend- 
ing on the (x, y) values) while in the large space the sums 
cluster between 21.2 and 21.4 for all of our (x,y) combi- 
nations. In the small space there are only seven states 
with J = 2+ T = 1 with the third state at about 18 MeV 
supplying 65-80% of the total isovector B(E2) strength. 

The great increase in the summed isovector B(E2) 
value when one goes to the large space is not due to the 
few lowest-lying J = 2 + T = 1 states. In the lowest part 
of table IB we see that in the ten lowest J = 2 + T = 1 



states in the large space we obtain an isovector B(E2)| 
sum that is 15-20% less than our sum in the small space 
and is less than 1/5 of the total summed B(E2) strength 
over the first 500 states in the large space, it will be in- 
teresting to see if there are one or a few very strong very 
high lying J = 2 + T = 1 states or whether this increase 
is a cumulative effect of many states. Another possibility 
is that we are getting into the region of spurious states. 

In the small space (and also for the ten low-lying states 
of the large space) the summed B(E2) increasing by 5% 
when the tensor force is decreased, increases by 22-33% 
when the spin-orbit is increased, and by 30-40% when 
both of the above effects take place. Very interestingly, 
however, in the large space, the summed isovector B(E2) 
is totally insensitive to the (x, y) values, varying by less 
than 1% among all the (x,y) combinations. 



IV. ADDITIONAL REMARKS 

The main problem addressed in this work in the fact 
that in a straightforward Op calculation with realistic 
G matrix elements obtained in a non-relativistic calcu- 
lation the computed value of B(Ml)f from ground the 
,/ = 1 + T = 1 state (0.89 /itjy) is much lower than the 
experimental value 2.63 This was already noted in 
the works of Zheng et al. j[5|,[l7). It should be empha- 
sized that in those works, as well as in this one, the single 
particle energies are calculated with the same interaction 
that is used between valence nucleons. 

In the small space the B(M1) can be enhanced by mak- 
ing the spin-orbit interaction stronger ||l5| , fi"7f and this can 
be justified by the Dirac phenomenology |L6| . The situa- 
tion with the tensor interaction is a bit more complicated. 
Since the Dirac phenomenology in its original version |l6| ] 
does not include pions, the tensor interaction is not af- 
fected. However, in the universal scaling mechanism of 
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G. E. Brown and M. Rho |ll| all mesons except the pion 
become less massive inside the nucleus as compared with 
free space. This leads to a weaker tensor interaction in- 
side the nucleus. This would be helpful for enhancing the 
value of B(Ml)t to the J = if T = 1 state. 

However, for other nuclei - 6 Li, A = 14 and 16 - 
another mechanism for causing the tensor interaction to 
become weaker was given by Zamick, Zheng and Fay- 
ache p0[ |. This was simply to allow AN = 2 or higher 
admixtures into the valence space, i.e. do larger shell 
model calculations or include the appropriate Goldstonc 
diagrams in perturbation theory. 

Whereas the quadrupole moment of the deuteron is 
positive, that of the J = 1+ T = ground state of 6 Li is 
negative (—0.082 e fm 2 ). In a (Op) 2 calculation without 
a tensor interaction, Q will be positive. With the real- 
istic Nijm II G matrix |^9| the value of Q is calculation 
to be —0.360 e fm 2 - too negative, indicating the tensor 
interaction is too strong. However, when 2hu> admix- 
tures enter Q becomes —0.251 e fm 2 , and with up to Ahuj 
admixtures Q is calculated to be —0.0085 e fm 2 . Thus 
it appears that high shell admixtures cause the tensor 
interaction to appear weaker in the valence space. 

In a similar vein, for the A = 14 beta decay 14 C [J = 
0+ T = 1] -> 14 N [J = 1+ T = 0] the experimental 
Gamow- Teller matrix element is very close to zero. It 
was shown by B. Jancovici and I. Talmi |5(J (see also 
D. R. Inglis |31j]) that one needs a tensor interaction in 
the (Op) shell space in order to get a vanishing matrix 
element. Again with the Nijm II interaction, which has 
a tensor part, the Ohui gives B(GT) = 3.967 - far from 
zero. With 2Tiui admixtures this gets reduce to 1.795. 
The effective tensor interaction in the valence space is 
weaker. If further we enhance the spin-orbit interaction 
by about 50% (as per discussion in the previous section) 
we can get B(GT) to be zero. 

In 16 O the isospin splitting of the J = + T = 1 and 
the J = Of T = states is 1.845 MeV. In the absence of 
a tensor interaction these two states are nearly degener- 
ate. With Nijm II the splitting becomes too large, 2.81 
MeV. In zero these two states are one particle, one hole 
state - admixture of the configurations (lSi/2 0p7y 2 ) and 

(0d 3 / 2 OP3/2 )■ When over and above these basic config- 
urations we allow 2tku excitations the energy splitting of 
T = 1 and T = is reduced from 2.81 MeV to 1.93 MeV, 
much closer to experiment. Again it appears that higher 
shell admixtures cause the tensor interaction to appear 
weaker. 

How do the higher shell admixtures affect the present 
problem 12 C? Logically, if we expect that higher shell ad- 
mixtures causes the tensor interaction to be weaker, we 
should do the large space calculation with y = 1, i.e. with 
the 'bare' tensor interaction. However, as seen in table |, 
the y = 1 calculations in the large space yields smaller 
values of B(M1) than do those in the small space. For 
example for x = the (small, large) results are 0.89 fi 2 N 
and 0.54 jj? N (recall that experiment is 2.63(8) For 



x = 1.5 i.e. an enhanced spin-orbit interaction the cor- 
responding values are 1.89 [M 2 N and 1.29 fi 2 N . 

The above results should not be totally unexpected. 
The isovector spin B(Ml) ff is proportional via an isospin 
rotation to B(GT), the Gamow- Teller transition from 
(N,Z)=(6,6) to (7,5) or (5,7), and for GT transitions 
it has been noted that a quenching factor is generally 
needed. A popular choice is 0.5 for B(GT) or 0.7 for the 
matrix element Mqt- Looking at B(M1) CT is table | we 
see that for x = 1, y = 1 it gets reduced from 0.58 fij^ 
to 0.31 n 2 N when 2Tiu) admixtures are allowed. This is 
almost the factor of two that is often quoted. 

So basically we are getting a quenching rather than 
the enhancement that we need. It would appear that 
the only way to get an enhancement towards the exper- 
imental value is to make the strength of the spin-orbit 
interaction x even larger and/or to invoke the univer- 
sal scaling argument of G. E. Brown and M. Rho pl| to 
make the tensor interaction weaker. Which, or how much 
of these two mechanisms should be invoked will require 
further investigation. 
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TABLE I. Comparison of the results of AN = (small space) and AN — + 2hui (large space) 
calculations of the value of B(M1) from the ground state of 12 C to the lowest (and most strongly 
excited) J = 1 + T = 1 state. The results are presented for different combinations of the spin-orbit 
(a;) and tensor (y) strengths in the realistic interaction of Eq. ([[J). 



small space 


X 


y 


E(MeV) 


B(Ml) a 


B(Ml) ; b 


B(M1) CT C 




1 


l 


13.60 


0.89 


0.03 


0.58 




1.5 


l 


13.08 


1.89 


0.04 


1.71 




1 


0.5 


13.33 


1.12 


0.01 


0.90 




1.5 


0.5 


13.11 


2.54 


0.00 


2.57 


large space 


X 


V 


E(MeV) 


B(Ml) a 


B(Ml); b 


B(M1) CT C 


AiV = + 27kj 


1 


1 


17.28 


0.54 


0.06 


0.31 




1.5 


1 


16.46 


1.29 


0.02 


0.98 




1 


0.5 


16.84 


0.88 


0.026 


0.60 




1.5 


0.5 


16.28 


1.96 


0.026 


1.83 


Experiment 11 






15.11 


2.63(8) 






a 9iTr — ^i9iv — 


0, g S n = 5.586, g sv = - 


-3.726. 








h g l7r = 0.5, g iv 


= -0.5, g sv 


= 0, g sv = 


0. 









c gu = 0, giu = 0, 9siz = 0.5, g sv = -0.5. 

d From Ref. ||. Also the isoscalar transition to the J = 1 + T = state at 12.71 MeV has a 
strength B(M1)=0.040(3) fj? N . 
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TABLE II. Same as table I, but for the summed strengths. In the small space, the sum is over 
all eight states. 



small space a 


X 


y 




B(M1) 




B(M1); 




B(M1) CT 


AiV = 


1 


l 


1.42 


0.60 


0.91 




1.5 


l 


2.29 


0.50 


2.07 




i 

i 


u . 


l.OO 


u.oy 


1 TO 

l.lo 




1.0 


n r * 
u . 


9 R^i 




9 89 
z.oz 


1 b 

large space 


X 


y 


B(M1) 


B(M1); 


B(M1) CT 


AiV = + 2nu) 


1 


i 


1.21 


0.56 


0.70 




1.5 


l 


1.77 


0.49 


1.43 




1 


0.5 


1.31 


0.56 


0.83 




1.5 


0.5 


2.27 


0.46 


2.085 


large space c 


X 


y 


B(M1) 


B(Ml)i 


B(M1) CT 


AiV = + 27kj 


1 


l 


1.61 


0.69 


0.97 




1.5 


l 


2.19 


0.62 


1.74 




1 


0.5 


1.57 


0.685 


0.98 




1.5 


0.5 


2.56 


0.57 


2.28 



a sum over all eight J = 1 + T = 1 states in the AiV — space. 
b sum over all first ten J = 1 + T = 1 states in the AiV = 2 space. 
c sum over the first 500 J = 1 + T = 1 states in the AiV = 2 space. 



TABLE III. Comparison of experimental B(E2) values to the first 2 + state in 12 C with shell model 
calculations for different combinations of spin-orbit (x) and tensor (y) interaction strengths. Effec- 



tive charges e v = 


1.5, e n — 


0.5 are 


used for AiV 


= calculations. 


Bare charges e ? 


, = 1.0, e„ = 0.0 


are used for AiV = 


= + 2huj calculations. 










X 


y 


E (MeV) 


B(E2) (e* fm 4 ) 


E (MeV) 


B(E2) (e A fm 4 ) 


Experiment* 1 






4.44 


39.9(2.2) 


16.11 


2.0(0.2) 


small space 


1 


1 


3.80 


80.6 


15.78 


2.6 


AiV = 


1.5 


1 


3.82 


75.3 


14.44 


3.1 


e p = 1.5, e„ = 0.5 


1 


0.5 


3.79 


80.1 


14.86 


1.2 












15.58 


1.6 




1.5 


0.5 


4.11 


72.9 


14.25 


3.8 


large space 


1 


1 


4.79 


23.2 


20.08 


1.6 


AiV = + 2Klu 


1.5 


1 


4.63 


21.7 


18.31 


1.1 


e p = 1, e„ = 










19.12 


1.2 




1 


0.5 


4.69 


23.2 


19.49 


1.5 




1.5 


0.5 


4.87 


21.2 


17.94 


2.3 



a From Rcf. W. 
b From Ref. fj. 
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TABLE IV. Comparison of the results in AN = (i.e. small) and AiV = 2 (i.e. large) spaces 
of shell model calculations for the sums of B(E2) transitions from the ground state of 12 C to 
J = 2 + T — states. The results are presented for different combinations of spin-orbit (x) and 
tensor (y) strengths in the realistic interaction Eq. (fil). 





X 


y 


isoscalar o(EjZ) b 


isovector ti\h,Z) 








J=U I = u — > J — I 1 — u 


J=U J=U — > J = Z 1 = 1 


small space 


1 


l 


82. 7 a 


3.39 b 


AN — 


1.5 


l 


78.4 


4.44 


total sum 


1 


0.5 


82.0 


3.54 


e p = 1.5, e„ = 0.5 


1.5 


0.5 


75.3 


4.80 


large space 


1 


1 


31.9 C 


2.91 d 


AiV = + 2hcu 


1.5 


1 


29.8 


3.55 


e p = 1, e„ = 


1 


0.5 


31.8 


3.08 




1.5 


0.5 


29.6 


3.87 


large space 


1 


1 


42. 5 C 


21.23' 


AiV = + 2hcu 


1.5 


1 


41.1 


21.40 


e p = l,e n = 


1 


0.5 


42.1 


21.24 




1.5 


0.5 


39.9 


21.36 



a sum over all eight J = 2 + T — states in the AiV = space. 

b sum over all eight J = 2 + T = 1 states in the AiV = space. 

c sum over all first ten J = 2 + T = states in the AiV = 2 space. 

d sum over all first ten J = 2 + T = 1 states in the AiV = 2 space. 

°sum over the first 500 J = 2 + T = states in the AiV = 2 space. 

f sum over the first 500 J = 2 + T — 1 states in the AiV = 2 space. 

g The experimental value to the 2^ state at 4.44 MeV is B(E2)=40 e 2 fm 4 (l3). 
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